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Abstract 

We extend the conformal gluing construction of |18| by establishing an analogous gluing 
result for field theories obtained by minimally coupling Einstein's gravitational theory with 
matter fields. We treat classical fields such as perfect fiuids and the Yang-Mills equations 
as well as the Einstein- Vlasov system, which is an important example coming from kinetic 
theory. In carrying out these extensions, we extend the conformal gluing technique to higher 
dimensions and codify it in such a way as to make more transparent where it can, and can 
not, be applied. In particular, we show exactly what criteria need to be met in order to 
apply the construction, in its present form, to any other non-vacuum field theory. 

1 Introduction 

One of the effective tools which have been recently developed and employed for the construction 
and study of initial data for solutions of Einstein's gravitational field equations is the method 
of gluing. This is a technique which has had a long and fruitful history in geometric analysis 
but which has only recently been successfully applied to general relativity. The idea of the 
gluing method (in its simplest form) is that, given a pair of sets of initial data which satisfy 
the Einstein constraint equations, we may use it to construct a new set of initial data which 
a) lives on a connected sum of the manifolds of the given sets of data, b) solves the constraint 
equations, and c) closely approximates the original sets of data on the parts of the new manifold 
which correspond to the original manifolds (i.e., away from the tubular "neck" of the connected 
sum). The work of Isenberg, Mazzeo, and Pollack ^S], shows that this sort of gluing can be 
carried out for fairly general sets of vacuum initial data. That work also describes a number of 
applications of gluing, such as producing multi-black hole initial data, adding wormholes to given 
sets of data, and showing that an arbitrary closed manifold with a point removed always admits 
both asymptotically Euclidean and asymptotically hyperbolic solutions of the vacuum constraint 
equations. 
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In this work, we show that the gluing results and the gluing applications which arc discussed 
for vacuum initial data in |18j can be extended to fairly general sets of non-vacuum data as 
well. We do this here for a number of special cases, including Einstein-Maxwell, Einstein- Yang- 
Mills, Einstein-fluids, and Einstein- Vlasov, as well as for any of these theories with a cosmological 
constant added. We also discuss the features which a general field theory should have if, when 
it is coupled to Einstein's equations, solutions of the corresponding constraint equations should 
allow gluing. 

As for vacuum data, the gluing procedure applied to non-vacuum data relies quite heavily on 
the conformal method for obtaining solutions of the constraint equations. Thus, after comment- 
ing in Section 2 on the general form of the constraint equations for non- vacuum field theories, 
we proceed in that section to describe the application of the conformal method to such the- 
ories. For the non-vacuum theories listed above-Einstein-Maxwell, etc-the conformal method 
leads to determined sets of (nonlinear) elliptic equations which, at least for constant mean cur- 
vature ("CMC") data, are readily analyzed for solubility. This is not true for all non- vacuum 
field theories; indeed, for any non-vacuum field theory which involves derivative coupling (e.g., 
the Einstein- vector-Klein-Gordon theory), the conformal method leads to equations which are 
intractable using known techniques. (See 1171.) Thus the gluing procedure which we use does not 
work for such theories. 

Note that in our discussion of the conformal method, we focus on the situation in which the 
initial data has constant mean curvature. We do this because, as with the vacuum case, even when 
the constraint solutions we are gluing together have non-constant mean curvature, the analysis we 
rely on to carry out the gluing is based primarily on the CMC version of the conformal treatment 
of the constraints. (See )19|.) 

While the details of the analysis differ from one non- vacuum field theory to another, the basic 
steps of the gluing procedure are largely the same for those non-vacuum field theories with the 
appropriate form for the constraint equations (in conformal form). Thus for these field theories, 
we can present a general discussion of these basic steps: (i) conformal blowup at the gluing points; 
(ii) connected sum of the manifolds and patching of the conformal metrics; (iii) patching of the 
non-gravitational fields, solution of the non-gravitational constraints, and deformation estimates; 
(iv) patching of the extrinsic curvatures, solution of the momentum constraint and deforma- 
tion estimates; (v) patching of the conformal factor, solution of the Hamiltonian constraint (in 
Lichnerowicz form) and deformation estimates; (vi) conformal recomposition of the initial data, 
forming the glued solution. We do this in Section 3. Also in that section, we outline the general 
analysis which leads to a proof that the gluing can be carried out for appropriate sets of initial 
data for these theories. The section culminates with a general statement of our gluing results for 
general theories. 

We discuss some of the details of gluing for various particular non-vacuum field theories in 
Section 4. Included are discussions of the field theories listed above. However, to avoid repetition, 
we focus primarily on two of them: Einstein-Perfect-Fluids and Einstein- Yang-Mills. 

Although the physically important versions of most of the field theories we discuss here are 
defined on 3-1-1 dimensional spacetimes, our results hold for arbitrary dimension. Hence we state 
most of our formulations and our results for n + I dimensional spacetimes where n > 3. 

A very different and important type of gluing construction has been developed by Corvino and 
Schoen ^U] , U21 and adapted and applied by Chrusciel and Delay [HI , [Z| ■ These results exploit the 
underdetermined nature of the constraint equations as opposed to using the conformal method 
to convert them into a determined system. This has lead to a number of remarkable results 
beginning with the existence of a large class of asymptotically Euclidean spacetimes which are 
exactly Schwarzchild near infinity. Recently, by combining these techniques with the results of 
[T^ and with the previous work of Bartnik 0, Chrusciel, Isenberg and Pollack [H], have 
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obtained a gluing construction of the type described here which is optimal in two distinct ways. 
First, it applies to generic initial data sets and the required (generically satisfied) hypotheses are 
geometrically and physically natural. Second, the construction is completely local in the sense 
that the initial data is left unaltered on the complement of arbitrarily small neighborhoods of the 
points about which the gluing takes place. Using this construction they have been able to establish 
the existence of cosmological, maximal globally hyperbolic, vacuum space-times with no constant 
mean curvature spacelike Cauchy surfaces. Except for the case of generic non-gravitational fields 
described entirely by an energy density function p and a current density vector field J (satisfying 
a strict energy condition p > \ J\), the Corvino-Schoen techniques have not yet been generalized 
away from the vacuum case. It is, however, expected that this can be done; it would then follow 
that the gluing theorems obtained for vacuum data in |5] would extend to non-vacuum data. 

We end this introduction by remarking that we have, for simplicity, restricted ourselves here 
to the consideration of initial data on compact manifolds (i.e., the cosmological setting). The 
gluing results presented here have analogous statements which are valid for either asymptotically 
Euclidean or asymptotically hyperbolic initial data sets. The required adaptations, which are 
similar to those discussed in detail in ^21; are left to the interested reader. 

2 The Constraint Equations and the Conformal Method 

We restrict our attention in this work to classical field theories which are obtained by minimally 
coupling a spacetime covariant field theory to Einstein's gravitational theory, and which have 
a well-posed Cauchy formulation. For such theories, if one is given a set of initial data which 
satisfies the constraint equations corresponding to that theory, one can always evolve to obtain a 
spacetime solution of the full PDE system. Our main interest here is primarily in the construction 
of solutions of the constraints. 

For the theories we are interested in here, the initial data consist of a choice of an n-dimensional 
manifold E" together with a Riemannian metric 7, a symmetric tensor K , and a collection of non 
gravitational fields which we collectively label JF, all specified on S". These non gravitational 
fields are usually, but not always, sections of a bundle over S". The constraint equations which 
these initial data must satisfy generally take the form 



where J is the current density of the non-gravitational fields, p is their energy density, and 
C denotes the set of additional constraints that come from the non gravitational part of the 
theory.^ Note that the first of these constraint equations is known as the momentum constraint, 
the second is often referred to as the Hamiltonian constraint, while the last are collectively labeled 
the non-gravitational constraints. 

As an example, for the Einstein-Maxwell theory in 3+1 dimensions, the non-gravitational 
fields consist of the electric and magnetic vector fields E and _B, we have p — }^{\E\'^ -t- |i3|^), 
J = {E X B)y, and we have the extra (non-gravitational) constraints div-yi5 = and div-y_B = 0. 

The system of constraint equations Q-®, vacuum or non- vacuum, is an underdetermined 
PDE system. The idea of the conformal method is to split the initial data fields into two sets of 
fields: the "conformal data", which is freely chosen, and the "determined data" which is to be 

^If we use T to denote the stress-energy tensor for the non gravitational fields and we use e±_ to denote the unit 
normal to the hypersurface E embedded in the spacetime solution generated from the initial data (S", 7, i^, JF), 
then J = —T{ej_, ) and p = T{ej_, ej_). Note that we have chosen units so that SvrG = 1 = c. 



C(^,7) 



divK -dtrK 



2p(^,7) 
0, 



(1) 

(2) 
(3) 
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found by solving the constraints. In the famihar vaeuum case [S], the conformal data consists of 
the manifold E", a Riemannian metric 7, a divergence-free trace- free symmetric tensor cr, and a 
function t, while the determined data consists of a positive definite function (p and a vector field 
W . With the conformal data (S", 7, cr, r) chosen, one determines ((/>, W) by solving the equations 

T? — 1 

div^ {VW) = — — <?!)«+2Vt (4) 

and 

- — + ^ |a + VW\l r"-"" - — rV^+' = 0, (5) 
q[n — 1) q(n — Ij ' qn 

where "DW is the conformal Killing operator, with coordinate representation 

VWab = VaWb + VbWa " 'lab^ cW\ (6) 

n 

and where q = -^^^ is a dimensional constant. If for a given choice of the conformal data one 
does determine W and (f> satisfying Q-©, then the fields 



7 = ^«7 (7) 
K = <l>-'^{a + VW) + -<j>'''y (8) 



n 



satisfy the vacuum constraints, consisting of Q-© with vanishing p and J. Note that the explicit 
form of the equations (@J and (O is determined by the form of the field decomposition, expressed 
in 0-©- The explicit choice of the form of the field decomposition is in turn determined to a 
large extent by the two identities (for 7 = ^^7) 

R^y -(p-'^-'^{q{n-l)A^^^ R^(l)) (9) 

(with q = being the imique exponent which avoids |V0p terms in ©) and 

\/^cl)-^Bab) = (10) 

which holds for any trace-free tensor B. 

The extent to which there exist unique solutions to equations and |SJ) for various classes of 
conformal data has been studied extensively; see [3] for a recent review. Here, we note primarily 
that while the issue is fairly well understand for constant mean curvature ("CMC") conformal 
data and near-CMC conformal data, very little is known more generally. Consequently, the 
earliest gluing results for the vacuum constraints ^H] pertain to CMC data sets, and later results 
rely primarily on CMC analysis. Hence we focus on CMC data sets here. 

A set of initial data has constant mean curvature if tr /iT = t is constant on E" . This condition 
significantly simplifies the analysis of the vacuum constraint equations because it decouples equa- 
tions and jSJ. Since = is a solution to © with vanishing right hand side, the analysis 
reduces to first finding a divergence-free trace-free symmetric tensor, and then solving (jSJ (often 
referred to as the "Lichnerowicz equation") with DW = 0. 

To extend the conformal method to the non-vacuum constraints, with non-gravitational fields 
T present, one needs to extend the field decomposition l|7Jl-lIHl) to JF. The chief criteria generally 
used to decide how to do this are the following ^7]: (I) In the CMC case, the constraint system 
should be semi-decoupled, in the sense that one can first solve the non-gravitational constraints 
independently of W and 4>, and then one can solve the momentum constraint for W independently 
of and then finally one solves the Lichnerowicz equation for (j). (II) The addition of the non- 
gravitational terms to the Hamiltonian constraint should not result in the Lichnerowicz equation 
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containing cither derivatives of (f> or those powers of (p which would lead to insurmountable 
difficulties in the subsequent analysis. 

To make these criteria more precise, let us presume that the field decomposition for the non- 
gravitational fields is defined by an action $ of the group of conformal factors {C^CS) under 
multiplication) on the set of matter fields, so that the physical fields JF, which together with 7 
and K must satisfy the constraints £^re given by = <I>(JF, (/)) where ^{T, 1) = J- and 

^{^{J^,(j)i),(j)2) = ^{J^,(j)i4>2)- Note that the data T is included in the set of conformal data 
(along with 7, a, and t); the explicit form of the action $ is to be chosen for each theory. 

In this language, the first of our criteria is satisfied so long as C, <&, and J satisfy the conditions 

C(<i>(.F,0),0^) = 0PC(^,7) (CI) 

and 

J(<i>(.F,0),0'7) = r'"V(^,7), (C2) 

for some number p. As a consequence of (|C1|I . if we choose the conformal data (7,(7, r,JF) so 

that C {T, 7) = 0, then whatever (p and W are determined to be, the constraint C {^T, 7^ = 

is satisfied by the physical initial data. As a consequence of (|C2|) . the conformal form of the 
momentum constraint is (compare Q) 

71 — 1 

div^ (VW) = (pt+'^VT + J{T, 7), (11) 

n 

which in the CMC case can be solved for W independent of 4>. 

Satisfaction of the second criteria depends upon the form of the term p($(J^, 0), ^'^7). It is 
crucial first of all that this quantity involves neither W nor any derivatives of (p. We formalize this 
by assuming that, for each choice of the matter field T and the metric 7, there exists a function 
njr^^ : R+ X S ^ M such that at any point p of S, 

<P'i+^2p{<^{T, 0), (/)'7) = n^Am.P)- (C3) 

Note that in practice, to ensure that this condition holds for a given field theory, it is generally 
sufficient that the non-gravitational fields do not involve derivative coupling. 

Assuming that condition (|C3|) holds, the Lichnerowicz equation takes the form 

</) - aiR-^(j) + ai |(7 + VW]^ + -a2T^(p'^"'^ + ai7i^^^(0) = 0, (12) 

where for convenience we set ai = , , ~ . J'^^.x and = — = -^^f-^, and also for convenience 
we suppress the dependence of njr .y on its second argument. To ensure that this equation is 
analytically tractable, it is important that for fixed 7 and T ^ nyr ,y(^(j)) is a monotonically decreasing 
function of (j). In practice, we find that this condition is often satisfied by nj^^^{4') being expressible 
as a sum of negative powers of 0, with non-negative coefficients. 

While the two criteria and their consequent conditions (|CH) - (|C3(I on $(^, (p) and its interaction 
with C, J, and p appear quite restrictive, it is shown in \\7\ that for most familiar physical fields, 
a choice of ^(JF, (p) which satisfies these criteria can be found. We shall discuss a number of such 
examples in this paper. Before proceeding, we illustrate how this works for a simple example: 
Einstein-Maxwell in three space dimensions . 

The conformal data for the Einstein-Maxwell theory consists of {Tr" ,^,a,T,B,E), where 
(I]",7, cr, r) are the usual vacuum conformal data, and both B and E are vector fields which 
are required to be divergence-free with respect to the metric 7. To satisfy the two criteria, we 
choose (^{B^.tp) = B''(p~^ and $(£"",(/)) = E°'<p-^. It follows then that (i) the extra constraints 
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div — and diY^E — are satisfied automatically so long as B and E are both divergence-free 
with respect to the metric 7; (ii) the momentum constraint takes the form 

div^iVW) = '^<P'^S/t + {E X B), (13) 

which is independent of (f) in the CMC case; and (iii) the Lichncrowicz equation takes the form 

We note that the extra term in (|14|l involves </) with a negative power and a positive coefficient, 
much like the |cr + LW^j term. As discussed in existence and uniqueness of solutions of 
P3|l - (|14|l in the CMC case then closely follows the pattern of the vacuum case. 



3 The Gluing Construction for Non- Vacuum Solutions 
3.1 Overview 

Before generalizing the gluing construction from |18| to apply to non- vacuum fields, we briefly 
summarize the original technique, modified to arbitrary spatial dimension n > 3. We start with 
an n-manifold E and a CMC solution (7, K) of the vacuum Einstein constraint equations (so 
K ~ a + ^7, with (T divergence-free and trace-free, and with t constant). We fix two points pi 
and p2 of E and a small radius R. Let Bj = Bii{pj) be the balls on which we will do surgery, let 
E* = E\{pi,p2} and let E* = Y\{Br{pi) U Br{p2))- The construction then proceeds as follows. 

• We first construct a conformally related metric on E* agreeing with 7 away from the surgery 
site and having two asymptotically cylindrical ends at the puncture locations. Let be 
a conformal factor equal to 1 on E2^ and equal to r^J'^ on Bj, where rj is the geodesic 
distance from pj, and where q = ^^^2 before. Then 7c = '(p'^'^l ^^e desired metric. 
Setting Oc = ipto' and Kc ~ ac + ^7 we see that (7c, K^) satisfies the momentum constraint 
and that satisfies the Lichnerowicz equation with respect to (7c,-?i'c)- 

• We next perform surgery on the cylindrical ends by identifying finite segments of length T 
to construct a family of topologically identical manifolds E^. To do this, we first construct 
maps from Bj\pj to the half cylinder (0,cx)) x by sending points at the ball radius 
rj to the cylinder length tj ~ — log rj + log R and by using Riemann normal coordinates 
to determine the projections onto S"^^. We then identify the finite segments {{tj,9) : < 
tj < T} via the map (ti, 6) 1-^ {T — ti, —9), resulting in the smooth manifold Et. Letting 
s = ti — T/2 = T/2 — t2, we denote by Qi.o the cylindrical segment {{s,9) : a — I < s < a + l} 
of length 21 centered at a. We use the shorthand notation Qi = Qi^ for centered segments, 
Q = Qi for a short collar in the middle, and Ct = Qt/2 foi' the entire cylindrical region. 

• We now construct approximate solutions (-^t , Kt , tpr) of the momentum constraint and 
the Lichnerowicz equation on Et by using the conformally modified solution away from 
the surgery site and using cutoff functions to piece together an approximation along the 
identified cylindrical segment. 

• In preparation for using the CMC-conformal technique, to map the approximate solutions 
to full solutions, we perturb the trace free part of Kt to obtain a constant trace second 
fundamental form Kt = o^t + ^7t which, together with 77^, satisfies the momentum con- 
straint. 
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• Finally, wc solve the Lichncrowicz equation using a contraction-map argument to arrive at 
a T parameterized set of solutions of the constraints which, for large T, is "close" to the 
original one away from the surgery site. 

The last step is quite delicate; the contraction-map argument demands some analytic condi- 
tions on the linearization Ct about of the Lichncrowicz operator A/r on Et- In particular, 
Ct must be surjective, and there must exist bounds uniform in T for the norm of its inverse on 
certain weighted Holder spaces. The argument also imposes some stringent conditions on the size 
of the error terms which arise from the earlier approximations and corrections. Substantial work 
in |18| is devoted to constructing and obtaining precise estimates for the perturbation of gt to 
gt that enter into this analysis. 

To extend the technique to include matter fields, we start as before with a CMC solution 
(7,/C, of the Einstein-matter constraints on an n-manifold S and a fibre bundle E over S. 
We also assume that we have chosen a conformal group action $ for the non gravitational fields 
which satisfies the criteria discussed in Section [2 

The topological step that constructs St must generally be supplemented with a construction 
of appropriate fibre bundles Et over St. Since each of the balls Bi is contractiblc, there exists 
a local trivialization over each of them. We can then use these trivializations to identify fibres: 
if qi G Bi is identified with (72 G B2 in the connected sum to create St, we can identify the fibre 
over qi with the fibre over q2 via our pair of fixed local trivializations. On the other hand, if 
q £ St is outside of the neck, then we take the fibre to be the one over q in E. The result is a 
smooth fibre bundle Ex over St- 

Together with the conformally modified gravitational fields jc and Kc, we define on S* the 
conformally modified non-gravitational fields 

Then as a consequence of the choice of $, we verify that (7c, Kc,J-c) satisfy the non-gravitational 
and momentum constraints on S*, and in addition ipc satisfies the Lichncrowicz equation corre- 
sponding to this data: 

A^, Vc - aiR^^i^c + aikc + - a2TV?+' + ain^,,^,(V'c) = 0, (15) 

To construct from {^c, Kc, J-'c) on S* a parameterized set of conformal data ['-^t , Kt , J^t) on 
St, we use a cut-off function procedure as in the vacuum case. That is, we first set (7T, Kt,J-^t) 
= {-fc,Kc,!Fc) on St\Q. Then on Q (recalling the definition of s in terms of T) we let x{s) 
be a cut-off function on M equal to for s > 1 and equal to 1 for s < — 1, and we define (i) 
7t = x(s)7i + (1 - X('S))72; (ii) OT = x{s)cri + (1 - x(s))cr2 in Q and Kt ^ (tt + ^7t (recall 
that T is a constant); and finally (iii) Tt = x{s)^i + (1 ^ x{^))^2- Note that ^i,(Ji, and !Fi (for 
i £ {1,2}) are all defined by identifying Ct with a subset first of Bi and then of B2. If we do 
the same with ipc, and then set ipr = x(*2 — l)'0i + x(^i ~ 1)V'2 in Ct and 1 outside, then we can 
verify (as discussed in further detail below) that for each value of T, (7T, Kt,^t) together with 
'0T constitute an approximate solution of the constraints, including the Lichncrowicz equation. 
Note that ipr = "01 + "02 on most of Ct- 

To go from these approximate solutions to a parameterized set of exact solutions, we pro- 
ceed as follows: First we perturb J-t, and thereby obtain a set of matter fields Tt which closely 
approximate J^c outside the neck, and which together with 7t solve the non gravitational con- 
straints globally on St. Then, we perturb Kt in order to obtain a set of symmetric tensors Kt 
which closely approximate K outside the neck and which together with 7t and J-t satisfy the 
momentum constraint globally on St. Note that both of these two perturbations involve solving 
partial differential equations for the perturbation terms (in order to satisfy the non gravitational 
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constraints and the momentum constraint respectively). The solvabihty of these equations (for 
the perturbation terms, given the approximate solutions) is an issue that must be addressed at 
this stage. Finally, provided that various error terms in the Lichnerowicz equation have been 
kept under control after substituting in 7^, Kt, J't, and V't, we use a contraction mapping argu- 
ment to show that there exists a solution ij^x to this equation with the conformal data (7^, 
Tt)', and we show further that away from the gluing region Q, the solution data (7T = i^xlT, 
Kt = ip^'^d'T + ^'i'T^TT, J-T ~ ^{^T,ijT)) approaches the original data {■-f,K,J^) arbitrarily 
closely. 

We discuss some of the details of these steps for generic non gravitational fields in the rest of 
this section, and discuss them for particular examples in section 4. 

3.2 Satisfying the Non-Gravitational Constraints 

Since we presumably have chosen the action of the conformal map <I> on the matter fields so that 
the non-gravitational constraints decouple from the others in the conformal representation, the 
first step in going from (jtjKtj^t) to (7T, Kt, ^t) is to choose !Ft so that C{!Ft,Jt) = 0. In 
obtaining J-t, we want it to be arbitrarily close (for T sufficiently large) to J-t on S^, in order 
to minimize the errors that are introduced into the other constraints. In particular, it is crucial 
for the gluing procedure that the error estimates HMip - l|M2p . HEip - l|E2() . and (|Nl(l - (|N4p described 
below be satisfied. 

The details of the construction of J-'t are field specific; we discuss a number of cases in section 
4. Here, we note what happens in the Einstein-Maxwell case: While the original E and B fields 
are divergence-free with respect to 7, and consequently the conformally mapped fields Ec and Be 
are divergence- free with respect to 7c, the fields Et and Bt constructed using cut-off functions 
are not divergence- free with respect to 7t (or any metric). Effectively the non-gravitational 
constraints are equivalent to this divergence-free property. We obtain new fields Et and Bt 
which satisfy the conditions div^^ Et — and div-^^ Bt = by carrying through the standard 
linear procedure. That is, we solve the linear equation 

A^y/iT = div-yy Et (16) 

for the scalar fiT, and then set 

Et^Et- VfiT; (17) 

the divergence- free condition immediately follows. We carry out a similar procedure to obtain Bt- 
Noting that the supports of div-y^, Et and of div-^^ Bt are contained in Q, one can carry through 
the analysis which verifies the estimates (|M1|I - (|M2(I . ljEl|l - ljE2|l . and ljNl|l - ljN4|l as detailed below. 

3.3 Repairing the Momentum Constraint 

With J-T determined, we next need to find Kt for which the momentum constraint is satisfied. 
We may do this by finding a symmetric trace-free (0,2) tensor i>t which satisfies 

div^y (;>t) = J(.?^t,7t) - div^y (TT- (18) 

If we can obtain a tensor j)^ which satisfies this condition, and if we then set Kt = o't + £'t + ^7t, 
then we do have a solution of the momentum constraint. 

To solve this equation for vt, we let V be the (77 compatible) conformal Killing operator 
on vector fields X, so VX = CxIt — ^div^^ {X)^t where Cx is the Lie derivative. Its formal 
adjoint V* is —div^j. and we set L ~ T)*!). If Wt is a vector field solving 

LWt = J{J^T, It) ~ div^y (Tt 
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(where we freely identify, via 77^, vectors and covcctors), then we can set Ot — —VWt to obtain 
a sohition to H18|) . and consequently a solution Kt to the momentum constraint. 

Besides obtaining i)T, we need to establish control of its size. Lemma ^ below, proved in 
for 3-manifolds S, provides this estimate in terms of the following Holder norm. 

Definition 1 Let \ \X\\k^aSi denote the Holder norm (computed with respect to the metric "fx) of 
the vector field X on an open subset Q ofS. Then we define 

ll-'^IU.Q = ||-''^|U.Q,S* „ + sup l|-'^IU.Q,Ql a- 

-T+l<a<i-l 

Lemma 1 Suppose there are no conformal Killing fields that vanish at the points pj ofT,. Then 
for T sufficiently large and for each X £ C'^'"(Et) there is a unique solution W G C'"'^^'"(S]t) to 
LW = X . Moreover, there exists a constant C independent ofW and T such that 

\\W\\k+2.a < CT^\\X\\k,^. 



The proof of Lemma ^ in general dimensions, which we skip now for the sake of exposition, is 
presented in Section 5. 

Since we wish i>t to be small, we therefore require that J{Tt, 7T)~div-j,y gt be small. Outside 
of Q, we have div^^ cry = J{!Fc,"lc)- Inside Q it is relatively easy to see that gt has norm and 
derivatives comparable to V'j'^'^ ^ e~""^/^. This motivates the following conditions on the matter 
field jFt- 

Definition 2 We say that Tt satisfies the momentum error estimates if for each k and a 
there exist constants C > and k > ^^^^ independent of T such that 

\\J{Tt,1t) - ^(•^c,7c)||,,„,s.\Q < Ce-^^ (Ml) 
and such that ( recalling that Qr is the centered collar of length 2r ) 

||J(^T,7T)|k-,o,Q. <Ce-'=^. (M2) 

If Tt satisfies the momentum error estimates it follows that 

||J(.FT,7T)-div^^ arlU^a < Ce""^. 

The threshold n > ^^^^ is the lower limit that will allow the subsequent analysis of the Lichnerow- 
icz equation to go forward. From the momentum error estimates and Lemma ^ we immediately 
obtain 

Proposition 1 Suppose Tt satisfies the momentum error estimates. Then there exists a tensor 
vt on TiT such that 

div^^ {vt) = J{Tt,1t) - Avv-iT (o't)- 

In particular, Kt = or + j^t + -^Tt together with jt o-nd Tt is a CMC solution of the momentum 
and non- gravitational constraints. Moreover, defining gt = vt + <^t, we find that there exist 
constants C > and k > ^^^^ independent of T such that 

||<TT-aT|U,a <Ce-'*^. (19) 



9 



The method of proof has been outhncd above; the details can be found m ifSj. We note that 
the existence of j>T-and therefore the existence of (jT-foUows even without Tt satisfying the 
momentum error estimates. The key point is that if we impose these estimates, the solution 
necessarily satisfies (|T^. 

3.4 Repairing the Energy Constraint 

Up to this point, we have constructed a CMC solution (jttKtjJ't) of the momentum and non- 
gravitational constraints; we also have an approximate solution ipT of the Lichncrowicz equation 
in terms of the conformal data (7T, Kt,J't)- Our goal is to find a perturbation rjT of ipT so that 
ipT + TjT solves the Lichncrowicz equation. 

Let Afr be the Lichncrowicz operator with respect to 77^, ax and J-'t on E^; we write 

Mt{iP) = A^^ij- aiRy^^j + ai |aT|'^ V"'"^ - a2TV+' + aifirii'), (20) 

where we abbreviate n^^ ^y(^) as hxi'tp). Similarly, we let Af denote the Lichncrowicz operator 
with respect to 7, a and on S. 

Since we will be working with the linearization of At, it is convenient to denote by ni^ the 
derivative of fix with respect to its last argument. Then the linearization Ct of A/t about tpx is 

Ct = \t - aiR-iT + o.i{-q - 3) |(7t|^^ V't''"'^ ^ "aCg + 1)t^V't + "^i "tI'/'t), (21) 

We will similarly denote by £ the linearization of M about 1. 

The contraction mapping argument we use here for the existence of a solution to the Lich- 
ncrowicz equation with data {-jtjKtj^t) requires that the error £t '■= Afri'ipT) decay faster as 
a function of T than a prescribed exponential rate; it also requires control of certain mapping 
properties of Ct- The proofs of these properties follow closely those given in taking into 
account the new features of A/y which result from the introduction of the term nx, as well as the 
use of higher dimensions. We now outline this argument. 

3.4.1 A Bound for the Error £t 

Corresponding to the momentum error estimates of Definition|21 we make the following definition. 

Definition 3 Recalling thatnTi'>pT) ■= ^^(jPt) depends implicitly on J-'t, and setting ndipc) ■= 

nj^^^j^{ipc) , we say that Tt satisfies the energy error estimates if for each k and a there exist 
constants C > and X > 1/q independent of T such that 

WM^t) - nMc)\l,c.,^^\Q < Ce-^^ (El) 

and such that 

\\nT{M\ko.,Q,<Ce-^^. (E2) 

Proposition 2 Suppose Tt satisfies the momentum and energy error estimates. Then for every 
k and a there exists a constant C > and a constant p > - not depending on T such that 

WEtWkc < Ce-P^. 

Proof: The proof follows the approach of the corresponding result in |18| . However, the compu- 
tations are more involved because we wish to prove the result for n > 3 and because of the effects 
of J-'t in 120|) . both directly via fiT and indirectly via (7^; hence we outline the proof now. 
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Wc divide Y,t into three regions, St\Ct, Ct\Q and Q, and estimate the various terms oi £t 
in each. For simpHcity of presentation, we only discuss estimates of the C" norm of £t; because 
of the exponential decay rates of the terms involved, the estimates for the higher derivative norms 
are essentially the same. It is also convenient to pick in advance the exponent 

• fx 1 1 2 ^ 3 

p := mm <A, — I — ,— ,2k — 1 

I 2 q q q 

One can readily verify that p > K In fact, the restriction 2k — 1 — | > i is the source of the 
choice K > 2^ in Definition |21 above. 

We first verify the bound in "EtXCt, where we have ipr = It = Ic and ctt ~ (^c- Since 
solves the Lichnerowicz equation with respect to {-jc, Kc, J-'c), it follows that in this region. 



where the definition of ric is analogous to that of ut- From Proposition ^ we know that 
|o't|^ — {o'tI^ < Ce^""^ in Et\Ct. From the energy error estimate (|F1|I it easily follows that in 
this region 

\£t\ < Cmax(e-''^,e-^^) 
< Ce-"^. 

Turning now to bounds in the region Ct, we first note the following easy estimates which hold 
on all of Ct- 

= ds^ + h + o(e-'^/^cosh{s)j (22a) 

< Ce-^/«cosh(2s/g) (22b) 
krl < Ce-"^/2(pogi^(2s/g))9+2 (22c) 
lar-crrl < Ce"'*^, (22d) 

where h denotes the round metric on the sphere. In the sub-region Q, we have (recall that 
^pT = i-'i + ip2 in most of Qt) 

£t = A/'t(^t) = (At - aiRT)ipi + (At - aii?T)V'2 + oi Io'tI^ V't^"^ ^ q-ir'^i'P^^ + aifiriipT) 
From (|22a|) and the definition of 7t it follows that At — aiRx = Ai — aiRi +0 (e^-^/^). Hence 

(At - aiRr) i'l = -a, ^Ar''' + «2t2^'?+' - a,n,{i^,) + O (e-(^+i)^) , 

where ni{ipi) := p($(J^i, i/'i), V?7i)V'?^^- Since ni{ipi) = ipf'^^ no{l) , for 7io(l) = p{J^-,l), we 
have the easy bound |ni(^/'i)| < Ce^^^^?)^. Also, since \ai\ < Ce"^-^ and tpj < Ce~"^^ it 
follows that 

\{At - aiRr) < C max (^e^^^^^^^ , e'^^^ 

an analogous estimate holds for '02. 

From l|22c|l and H22d|l we determine that 

|aT|'<Cmax(e""^,e"2'^^), 
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and hence from (I22bll we have 



From (|22H|) we also know that r^ip^^^ < Ce"(^+i)^ and from that |nT('(/'T)| < Ce~^^. 
Hence we verify that the estimate 

\£t\ < Ce-P^ 

holds in the region Q. 

The remaining region Ct\Q has two components, ■* — [— r/2, 1] x 5"^^ and c'^^ = 
[l,T/2] X 5"^^. By symmetry it suffices to prove the bound on just one component. In c'^'' we 
have (recall that V't = "02 + XiV-"! in C't ) 

= J^tH^t) = (At - aiRT){xii'i) + a-i {\<7t\'^ V'^'^"^ - l^rl^ V't"^) 

-a2T2 (V^fi - 0|,+i) + fli (nT(V'T) - ric(V'c)) • (23) 
The last two terms of H23|l are easy to estimate. From ()E1|I we know 

\nT{M-n,{i^,)\<Ce-^^. (24) 

Since tpr = + Xii^i ^'^d since xiV'i = ( e^"^" ) we have 



< Cmax(^e-f^,e~(^+?)'^ 



(25) 



We now turn to the second term of the right hand side of if^ . From (|22bp . H22c|) and (|22dp 
we have 

IcttI' - ^tI'I 0?'"' < Cmax (e-(''-i)^, g-C^'^-i-f )^) . (26) 

Also, since 



we have from (|22cf) that 



I e 9 , e ^ I 



-pT 



(27) 



(28) 



From (EHJ and it then follows that 
Hence it remains to estimate 

(At - aiRT){xii>i)- 

Since xi = 1 except near s ~ T/2 and since ^pl ~ e^T/q-2s/q it follows that 

Xi^i^e-^/''-'^/'^ + 0[e--^^). (29) 

Letting 70 be the round metric on the cylinder it follows from H22a|) and from calculation of the 
scalar curvature of a round metric that 



At - ciRt = Ao - ai(n - l)(7i - 2) + O fe'^^/^ 



(30) 
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Since (Ao -ai{n- l)(n - 2)) e-^/9-2^/« = 0, we see from ^ and (EOJ) that 

(At - aii?T) (xiV-i) <Cmax(e-(^+^)^,e-'^'^^ . (31) 
Hence from if^ . (j2Hll, and (jSU it follows that in ^ (as in the rest of St) 

□ 

3.4.2 Mapping properties of the linearization of Afr 

The goal of this section is to verify certain properties of the linearization operator Cx- specifically, 
that it is invcrtible and that the norm of its inverse on certain weighted spaces is bounded 
independent of T. These properties, which hold in the vacuum case provided K ^ 0, depend on 
the form of the matter term n{'^^!F, ip) and its derivative n' with respect to ip. To guarantee these 
properties, we add the following assumptions to the momentum error estimates and the energy 
error estimates already made: 



• For each p G S, 

• As T — ^ oo, 

• As T — > oo, 



n^,^(l,p)-n^_^(l,p)>0. (Nl) 

ll"T(V'T)-<(V'c)||fe,„,S,\Q^O. (N2) 
II^T(V'T)IU.„n, -0. (N3) 



We use property (|N1|) to establish that C (the linearization of A/" on S about -0 = 1) is invertible, 
and we use properties l|N2|l and l|N3|l ensure the convergence of Ct to known operators away from 
the middle of the neck and on the middle of the neck respectively. 

To see that C is is invertible, we note that 

£ = A^ - ai (r^ + (g + 3) + ^{q + 1)t^ - n'^Jl) 

Since (7, cr, t, J^) is a solution of the constraints, we can write in terms of 7, a, r and J- to 
obtain 

£ = A^ - ai (^{q + 4) \a\l + g^r^ + („^,^(1) - n^,^(l))) . (32) 



As a consequence of our assumption (|Nip . we know that the term in this expression which involves 
is non-negative. Moreover, if the term in parentheses is not identically zero then it follows from 
the maximum principle that C has a trivial kernel on C'^'+^'"(ri). This is true if cither K ^ 
or — n'jr ^{1) ^ 0; we will henceforth assume this non-degeneracy condition. Note that 

it often holds that — n^^(l) > 0. In this case, since n;F_^(l) — 2p(jF, 7), the non-degeneracy 
condition holds if either A' ^ or p{T, 7) ^ 0. 

We now show that for sufficiently large T, Ct also has trivial kernel. We also want to control 
the norm of its inverse, which we can do in terms of the following weighted Holder space norm: 
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Definition 4 Let wt be an everywhere positive smooth function on Y,t which equals e~"^l'^ cosh(^) 
on Ct, which is uniformly bounded away from zero on and which is equal to one on ^2;^. 
For any 5 G K, and any (p G C'^'"{'St), we set 

ll0l|fc,Q,5 = \\wT^(l>\\k,a, 

and we fei C^'"(I]t) denote the corresponding normed space. 
Proposition 3 Fix any 5 G R. For T sufficiently large, the mapping 

is an isomorphism. 

Proof: Since the spaces are identical as sets, and since the associ- 

ated norms arc equivalent to each other for fixed T and (5, it is enough to show that Ct is an 
isomorphism for (5 = 0. Since Ct is Fredholm of index zero on C'^+^'"(S]t), we need only show 
that, for T sufficiently large, it has trivial kernel. 

Suppose not, so that there exists a sequence of increasing parameters Tk and functions rjk (not 
identically zero) such that CT^Vk = 0. We assume without loss of generality that sup \ i]k\ = 1- 

Suppose now that supl?]^! is bounded away from zero on S* for some r < R. As a conse- 
quence of elliptic regularity, together with property ljN2p . we conclude that there exists a nonzero 
function 77 on S* such that Ccrj = 0. Using the conformal covariancc property A-y^ r/ — aiR^^rj ~ 
^P'^ (A-y {•q/^c) — aiRjr]/ipc), we conclude that r]/'iljc satisfies 

A^ - aiR^ -{q + 3)ai |c7|^ + ^2(9 + 1^ - ai^-'^n'M] (v/i'c) = 0. 

Now, it follows from the definition of n that nc{ip) = ip^'^^n{ipip^-^), and hence that ip'^''n'^{ipc) = 
n'(l). So rj/tpc solves C{ri/tpc) = on S*. Since 77 is bounded on S* and since ipc decays like 
rj" ^■'''^ near pj, it follows that i]/ipc is less singular than rj" at pj and hence extends to a 
nontrivial solution of C{i]/'ipc) = on all of S, which is a contradiction. 

If on the other hand rjT converges uniformly to on S* for any i? > r > 0, we can consider a 
sequence of functions on increasing finite length sections of the cylinder E x 5*""^ by translating 
r]T so that its maximum occurs at s = 0. Then it follows from elliptic regularity and the properties 
()JN2|) and 1N3() that we can extract a subsequence that converges in on compact sets of the 
cylinder to a nontrivial bounded solution of the equation 

Since there are no such solutions, we again have a contradiction. □ 

Let Ot denote the inverse of Ct on Cg^'^'" , which exists for T sufficiently large, as a conse- 
quence of Proposition 3. The proof that Qt has bounded norm as T goes to 00 is, using properties 
HN2|I and (|N3p . identical to the one appearing in ^Sj, and hence we will not repeat it. 

Proposition 4 IfO<S<l, then the norm of the operator Qt *s uniformly bounded as T ^ 00. 
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3.4.3 Deforming the conformal factor if^T to a solution 

We deform ipx to a true solution ipT = V't + Vt of the Lichnerowicz equation using a contraction 
map. That is, we wish to solve 

XriipT + Vt) = 0. 

To do this we introduce the quadratically vanishing nonlincarity Qt defined by 

AfriipT + v) = J^tH't) + Ct{v) + Qriv)- 
Then ipx + ryr is a solution of the Lichnerowicz equation if and only if rjT satisfies 

Vt = -Qt (St + QTivT)) , 
and hence is a fixed point of the map Tt defined by 

rriv) = -Gt (St + Qt{v)) ■ 

The key estimate which leads to a verification that 7^ is a contraction map is that, for r) near 
zero, the nonlinear operator At differs from its linear approximation Ct by a quadratically small 
amount. That is. we wish to show that there is a constant C independent of T such that 

\\QTmk,a,S<CM\l^^,. (33) 

For the terms in Qt involving a and t, this estimate follows from the explicit form of their 
dependence on 77. To argue this, let us for convenience define 

wt{x) ai\aT? {x-^-^ - aar^ai^+i) . 

We then have 

Qt{ii) wriipT + 1]) - ujTiipr) - w'rp{i}jT)'q + riTii^T + v) - nrii^r) ~ »^t(V't)'?- 
The terms involving wt can be rewritten in integral form as 

wrii'T + Vt) - wxiil^T) - w'rp{ipT)riT = / (vt - s)w'^(ipT + s) ds; 

JQ 

hence an estimate of the type (|33ll for these terms readily follows from appropriate controls on 
Wxi'ipT + s). The only serious difficulty stems from the negative exponents appearing in the 
definition of wt- Following |18j . we require that there is a constant c < 1 independent of T such 
that I77I < cipT- This restriction ensures that 1 + v/'i'T remains uniformly bounded away from 
zero and a straightforward computation using this fact shows that an estimate of the form H33|) 
holds for the wt terms. 

Since there is no a priori form for the dependence of the matter terms ut on ^y, to complete 
the estimate for Qt we need to make a further assumption on tit (thus extending our list of 
assumptions from section TS. 4.2(1 : 

• There exist constants C > and c < 1 independent of T such that 

WnriipT + v)- nriM - "tIV't)??!!^,^,^ ^ (^M\l,a^5 (^4) 
whenever \r]\ < c\iPt\- 
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With the estimate conditionaUy in hand, we now seek an open neighborhood about 
on which Tt (as defined above) is a contraction map. This neighborhood must be mapped to 
itself under the action of Tt and it must be smaU enough so that (for all T) it only contains 
functions rj such that \'q\ < cipT and hence the quadratic estimate holds for it. Fixing i5 > 
we find from Propositions [21 and 01 that there are constants C and M respectively such that 
||t;(£T)||fe a 6 ^ CA/e("^+9)^. It is therefore convenient to set 

M, = {ueCl''' : ||w|U,a><i^e(-'+f)^}, 

and we note for future reference that \i v — 2CM then Gt{£t) G B^. 
Now for J] £ Ml, we see that inside the neck 

\ri\ < i^e("^+f)^VT~VT 

whereas outside the neck 

Since (we recall from the energy error estimate) A > -, it follows that for any fixed values of 
c < 1 and v there exists a T large enough so that if 77 € Bj^, then \r]\ < cipT in all Sy. Combining 
this inequality with estimate we find that there exists large enough T so that 

for 771 and 772 in B,y. 

Setting ~ 2CM we have G{£t) S B^^- Moreover, setting 771 = 7/ and 7/2 = in H3.4.3|) it 
follows that Q{Qt{v)) & I^f for r] G M,, and for sufficiently large T. It follows that Tt maps 8,^ 
to itself for sufficiently large T. Equation H3.4.3|l also implies that Tt is a contraction map on Mi, 
for T sufficiently large. (See QHl for further details.) 

If is a contraction map, it follows that it has a unique fixed point rjT. It then follows that 
ipT = + Vt satisfies the Lichnerowicz equation 

Afriipr) = (34) 

with A/t defined (for the glued data ^t,o't,^t,t) as in (|20|l . Moreover, the difference between 
the approximate solution ipx and the exact one ipT is contained in the space B^, which we see 
from its definition gets arbitrarily small as T gets large. We are thus led to our main result: 

Theorem 1 Let (S", 7, K, T") he a smooth, constant mean curvature solution of the n-dimensional 
constraint equations for an Einstein-matter field theory which satisfies the two criteria of Section 
2. Let S" be compact, and let pi and p2 be a pair of points in . We assume that this solution 
data satisfies the following conditions: i ) The metric is nondegenerate with respect to pi and p2 in 
the sense of lltl^ ^ , and either the quantity K or the quantity nj^^^(l) — 77,^ ^(1) appearing in con- 
dition \N1\) is not identically zero.^ ii) The momentum error estimates (Definition 2), the energy 
error estimates (Definition 3), and the matter term estimates fNl\} - ![N4}^ all hold. Then, for T 

^That is that we assume there are no non-trivial conformal Kilhng fields on which vanish at either pi or p2. 
^Note that for most Einstein- matter field theories, it follows from the momentum constraint that if is not 
identically zero, then K is not identically zero. 
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sufficiently large, there is a one-parameter family of solutions {T,1^,Tt, Kt, Ft) of the Einstein- 
matter field constraints with the following properties: The manifolds ( all diffeomeorphic to 
each other) are all constructed by adding a handle or neck to E", connecting the two points pi and 
P2. On the region of Yi^ outside of the neck, the fields {Tt, Kt, Ft) approach arbitrarily closely 
to (E",7, JF) as T tends to infinity. 

If we compare the hypotheses of this theorem for gluing solutions of the Einstein-matter field 
constraints with the hypotheses of the corresponding gluing theorem for the Einstein vacuum data 
(See Theorem 1 in jl8)). it appears as if there are far more conditions which must be satisfied by the 
Einstein-matter data. We note, however, that for most of the Einstein-matter theories of physical 
interest (See Section 4), all of the conditions except nondegeneracy are satisfied automatically by 
solutions of the constraints. That is, the nature of the matter fields and how they couple into 
the Einstein-matter constraints for most such theories guarantees satisfaction of most of these 
hypotheses. We shall see this in the discussion of example theories in the next section. 

3.5 Gluing in the Presence of a Cosmological Constant 

What happens if we seek to glue solutions with a cosmological constant A present? Wc find that 
so long as the constant a2T^ — aiA is positive (or zero in certain appropriate cases), our results 
are unchanged; otherwise, gluing via the methods we describe here generally cannot be carried 
out. 

To see this, we recall that since the field equations with A present take the form G = T + Ag, 
it follows that the momentum constraint H4U|) is unchanged by the presence of A, while the 
Hamiltonian equation (|41|) is changed to the following: 

i?^-|A'|2 + (trii:)2 = 2p + 2A. (35) 

Then, since we require A to be invariant under the conformal action $ (otherwise we could not 
guarantee that it would remain a constant) wc find that the Lichncrowicz equation with A present 
takes the form 

Aa - aiR^cj) + ai\a + VW\l 0-^-3 _ ^^^^2 _ aiA)(f>''+^ + n^^^{c^) = 0. (36) 

The pairing of A with which we see in H36|l persists throughout the gluing analysis. Since it is 
crucial to this analysis that a2T^ — aiA > 0, our stated result follows. 

4 Applications to Example Einstein-Matter Field Theories 

The work done in Section 3, culminating in Theorem 1, provides criteria for determining whether 
gluing can be carried out for given solutions of the constraints of a specified Einstein-matter field 
theory. While these criteria involve the data of the particular solutions, they depend predom- 
inantly on the features of the specified theory. In this section, wc discuss some example field 
theories for which those solutions of the constraints which satisfy the mild nondegeneracy condi- 
tions stated in hypothesis (i) of Theorem 1 automatically satisfy hypothesis (ii) as well, and can 
therefore be glued. 

4.1 Einstein-Perfect Fluids 

The Einstein-perfect fluid field theories are especially simple from the point of view of gluing 
because they do not require any extra bundle structure, they introduce no additional constraints 
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beyond the Hamiltonian and the momentum constraints of the vacuum theory, and the fluid fields 
can be readily scaled in such a way that the conformal method analysis is without complication. 

The additional fields needed for an Einstein-perfect fiuid field theory are the fluid energy 
density C, the pressure p, and the fluid velocity u (which satisfies the condition g{u, u) = —1). In 
terms of these quantities, the stress-energy tensor for the field theory is given by 

T = Cu(»u+p{g + u(»u); (37) 

and the field equations consist of 

G = T (38) 

together with 

divT^O, (39) 

which is a consequence of (|37|l . The system is complete, with a well-posed Cauchy problem, once 
an equation of state specifying the pressure p as a function of the density, p = ^(C)i has has been 
given. 

Regardless of the choice of equation of state, initial data for a solution of the Einstein-perfect 
fiuid equations consists of the vacuum initial data (S, 7, K) together with a scalar field C and 
a spatial vector field v (the spatial projection of the fluid velocity u). The constraint equations 
take the form 

divK -dtrK = J = ~{l + \v\^)^{C + F{C))v (40) 
i?^-|X|2 + (trA')2= 2p =2C + 2(C + F(C))|t;p. (41) 

To extend the conformal method to these Einstein-perfect fluid constraints, we need to define 
the conformal action $ on the matter fields. While one could do this working directly with the 
fiuid initial data fields C and v, the analysis is much simpler if we work instead with the composite 
functions p and J. If we are to do this, we need the following result: 

Lemma 2 The map from {(,v) to {p,J) given by p{C,v) = C + iC + F{C))\v\'^ o.'i^d by J{C,v) = 
— (1 + |vp) 2 ((^ 4- i^(C))v is invertible, so long as we assume the physical conditions p^ > |Jp, 
P = F{C) > 0, F'iC) < 1, and C > 0. 

Proof: To prove this lemma, let us presume that we are given (p, J) , and that we want to solve 
for (C,'y)- We first notice that it follows from equation 1)40(1 that J and v are parallel. Hence it is 
useful to calculate |Jp = |vp(l -I- |up)(C + F(C))^ = {p - C){p + F{C)), and then focus on solving 
for (C, |up) in terms of (p, | Jp). 

Let us fix p and define the function Hp{() := {p — C,){p + F{Q). Note from the form of the 
equation p{C,v) = C + (C + F{())\v\'^ and the hypotheses that C ^ and F{() > be non- 
negative that if a solution exists, we must have < ^ < p. So we seek to solve the equation 
Hp{C) = |Jp (fixed J) for C with < ( < P- Using the conditions F(C) > and > |Jp 
we find Hp{0) = p{p + F{Q)) > p^ > | Jp and also Hp{p) = < | Jp. We conclude from the 
continuity of Hp that there exists ( S [O,/?] such that ffp(C) = Calculating then that, as 

a consequence of the condition F'{C) < 1, we have -ffp(C) < 0, it follows that there is a unique 
C(p, |Jp) which satisfies -ffp(C) ~ \ J\'^- We may use this result together with equation H41|) to 
solve for |wp(p, |^P), and then finally we may use equation (|40|l to solve for v{p, J). □ 

As a consequence of this Lemma, we may now treat p and J as the initial data variables for 
the fiuid field, in place of C and v. Hence, to extend the conformal method to the Einstein-fiuid 
field, we define conformal action <I> on p and J: we set $(p, </>) = p0" 29-2 ^^ j^ ^-j _ J0^«~2. 

This choice has three consequences: i) The quantity ^■i'''' is conformally invariant, and hence 
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the satisfaction of the dominant energy condition by the fluid field is also conformally invariant, 
ii) The quantity J satisfies the condition (jC2p . and hence the conformal momentum constraint 

n — 1 

div^ {LW) = 0«+2Vr + J (42) 

n 

is independent of (/) so long as t is constant, iii) The matter-dependent term in the Lichnerowicz 
equation 

Aa - aiR^(l) + ai\n + VW\l - asr^^^+i + aip0-9/2"i = 0, (43) 

like the \a+'DW\ term, has a positive sign and contains a negative power of (p; hence its role in the 
solvability analysis of H43(l is essentially the same as the \a -\-'DW\'^ term. Thus we find, with this 
choice of the conformal action on p and J, the conformal method works for the Einstein-perfect 
fluid constraints more or less the same as it works for the Einstein vacuum constraints. This is a 
prerequisite for carrying out the gluing of solutions of the constraints as described in Section 3. 

Let us say that we are given a constant mean curvature solution (E, 7, K, C, v) of the Einstein- 
perfect fluid constraints, together with pair of points pi,p2 S S at which we would like to carry 
out a gluing operation. We first rewrite the initial data in the (S, 7, if, p, J) form; we work 
entirely in this form until the end, at which time we may invert to recover the glued solution in 
(S,7,if,C,i') form. 

The steps which take us from (S, 7, if, p, J) to the preliminary glued data sets (7^, Pt, Jt) 
are straightforward, as described in Section 3.1. With no non-gravitational constraints C present, 
we set Pt = Pt and Jt = Jt, and we next proceed to repair the momentum constraint, as in 
Section 3.3. To continue with the gluing at this stage, we need to check the conformal Killing 
field ( "CK" ) nondegeneracy condition (described in footnote |21 above) and we need to verify 
the momentum error estimates (|M1|) - (|M2|) . The CK nondegeneracy condition is (as with the 
vacuum case) a mild restriction on the class of solutions of the constraints which admit gluing. 
The momentum error estimates, on the other hand, follow immediately from the choice of the 
conformal action $ on (p, J); they are automatic for any solution: We verify 



I^T- J)||,,„,5:,\Q=0 (44) 

and 

\\JT\\k.a,Q2 < Ce^^, (45) 

thereby confirming that these estimates hold. Consequently the estimate (|19f) holds as well. We 
emphasize that this is true for any chosen solution (E, 7, if, p, J) of the Einstein-perfect fluid 
constraints which satisfies the CK nondegeneracy condition relative to the chosen gluing points 
pi and p2 ■ 

We next need to check that the energy error estimates (|E1|) - (|E2|I are satisfied. Noting 
that, for the Einstein-perfect fluid theories, rir(7T,PT, Jt,iPt) = PtV'"''^^"^: these estimates are 
readily verified for any choice of data, much like the momentum error estimates just discussed. 

To ensure that the linearized Lichnerowicz operator Ct is invertible for the chosen data 
(S, 7, if , p, J), we need the quantity njr_^ to satisfy condition IjNip . This follows immediately 

from the expression nT{'^T,PT, Jt,'>1-'t) = pt^~'^^~^"^, regardless of the choice of data. If we 
are working with data on a closed manifold, we also need one or the other of the conditions if 7^ 
or p ^ to be satisfied. This is the only restriction beyond the CK nondegeneracy condition that 
we need to make on the data itself. Note that it is analogous to the restriction imposed on the 
data for gluing in the vacuum case. 

The remaining assumptions which need to be verified are ljN2p . IN3|) and ljN4p . The first 
two of these essentially follow from the conformal invariance of 2l4iA. The last of these IS a 
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straightforward calculation. We thus conclude that any set of constant mean curvature initial 
data satisfying the CK nondegeneracy condition relative to pi and p2 and having either K ov p 
nonvanishing can be glued at pi and p2. 

4.2 Einstein-Maxwell and Einstein- Yang-Mills 

We first briefly review the Yang-Mills field theory (on a fixed background spacetime) to establish 
our notation. Let be a vector bundle with compact structure group G over a Lorentz manifold 
(M, rf). The sub-bundle of End(A^) consisting of those transformations associated with the adjoint 
representation of G is denoted ad(A^); note that each fibre of ad(A^) is isomorphic to g, the Lie 
algebra of G. We fix a bi-invariant metric on G, which induces a metric on the fibres of ad(A1). 
If I? is a connection on Al, then its curvature Fx> = 2?^ is a 2-form taking values in ad(A^), and 
therefore is a section of Vl^{a.d{M)). 

A solution to the Yang-Mills equations is a connection T) such that T>*Ft> = 0, where V* = 
^_^^dim(M)+i ^ where * is the Hodge star operator. (Note that to define this operator 

we need the specified metric 77.) One readily verifies that the Yang-Mills system has a well-posed 
Cauchy problem, for which the initial data consists of a G vector bundle V over a Riemannian 
manifold S with a connection D and a section E of ri^(ad(V)). The variable E plays the part of 
the time derivative of D and necessarily satisfies the constraint equation D* = 0, where D* is 
defined analogously to V* . In the discussion below, the curvature of D appears; we denote this 
curvature by Bd , recalling its familiar role in Maxwell's equations. Note that Bd is a section of 

Maxwell's theory is a special case of Yang-Mills theory, in which the structure group is chosen 
to be U{1). Since g in this case is naturally isomorphic to M, E and Bd are real- valued differential 
forms. In three dimensions, E and *-Bd sue the covector fields corresponding to the electric and 
magnetic fields. 

Minimally coupling Yang-Mills to Einstein we obtain the Einstein- Yang-Mills theory (with 
Einstein-Maxwell as a special case), for which the field equations are 'D*Fxi = and G = T, 
where V is now the covariant derivative corresponding to the gravitational as well as the Yang- 
Mills connection, and where the Yang-Mills stress-energy tensor is (in component form) Tap = 
PaPtj.13 ~ ^gapF^'^ Fpiu- Using techniques similar to those discussed in one readily verifies 
that the Einstein- Yang-Mills field equations have a well-posed Cauchy formulation. A set of initial 
data for Einstein- Yang-Mills consists of a G- vector bundle V over an n-manifold S, with Einstein 
data 7 and K , together with Yang-Mills data D and E, all satisfying the coupled constraint 
equations 

diviC- dtriC = 2(-l)"+H (£;,*Bd) 
i?^- + (trA')2 = |i;|2 + |BD|2 

V>* E = 0. 

Here the bilinear map (•, •) on sections of f2*(ad(V)) is the one which is induced from the bi- 
invariant metric on g, and the norm on sections of 0'^(ad(V)) is defined by \W\'^ = * < W,*W >. 
For Maxwell's equations, the equation D* E' = is the familiar constraint div E = (the magnetic 
field constraint div B = is absent since it is automatically satisfied when Maxwell theory is 
formulated as a special case of Yang-Mills theory) . 

The choice of the conformal action $ on the Yang-Mills fields D and E is essentially determined 
by the criteria discussed in Section 2. In particular, noting that jc = ''/''t implies that the Hodge 
star operator acting on fc-forms transforms via *c , we find that in order to satisfy 

(jCl|) and thereby prevent the appearance of i/j in the conformal version of the non-gravitational 
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constraint D* E ~ 0, wc require that ^{D,^) = D and $(£',■0) = ip^^E. It then foUows that if 
the conformal data satisfy D* i? = 0, then the reconstituted data satisfy D*E = 0. For MaxweU's 
equations in three dimensions, the conformal action on the connection $(Z3, ^p) ^ D is equivalent 
to a conformal action on the magnetic covector field B given by V') ~ ip^'^B. 

With this choice of the conformal action, we may proceed to check that (|C2|) and (|C3p hold 
as well. We calculate 

= ^ (E.^Bu) , 

thus satisfying condition ljC2|l : and we calculate 

{mE,i,)\i + ibdI?) - + v^-'ii?Dp, 

from which (|C3|) follows, with 

To apply the gluing construction, we start with a solution {j, K, E,D) of the Einstein- Yang- 
Mills constraints on a G vector bundle V over E. From the induced vector bundle Vc over S* with 
connection Dc we construct the vector bundles Vt over as described in Section 3 by fixing 
local trivializations on the balls Bj and using them to identify fibres over identified points in the 
connected sum. We then construct a connection Dt on Vt given by Dc except over Ct where 
Dt = Di — x(s)) D2 (here x is the cutoff function used to define 7r). Hence Dt = Dc 
except over Q. In terms of the classical notation for Maxwell's equations in three dimensions, the 
construction of this connection results in setting the magnetic field to curl(x(s)Ai -I- (1 — x{^))^2) 
over Ct, where each Aj is a vector potential for the magnetic field over the ball Bj. 

To construct Et, we first note that we may define local sections Ec, Ei, and E2 of the bundle 
ad(VT) restricted to T,\Q, Ct, and Ct respectively by identifying ad(£T) with ad(£c) over the 
appropriate subdomains. We then set Et ~ Ec over S\(5 and we set Ex ~ XiEi +X2E2 over Ct, 
where xi = x(^2 ~ 1) and X2 = x(^i ~ !)■ Thus Et is pieced together in the same way that the 
conformal factor ipx is (see Section 3.1) and one has Et = Ei+ E2 over most of Ct- It would be 
more convenient to stitch together Et over Q alone. However, doing so results in unacceptably 
large error terms and would result in a failure to satisfy the momentum error estimates. 

The next step in the gluing construction is the repair of the non-gravitational constraint: we 
need to find a section Et of ri^(ad(VT)) which satisfies D^ Et = 0. To do this, we seek a section 
/iT of ad(VT) that satisfies 

Y)rp Dt fJ'T — — T)rp Et- 

Then Et = Et + Dt A*t is the section we need. The operator DJ, Dt is self-adjoint and elliptic; 
therefore since — D^ Et is orthogonal to the kernel of D^ Dt, we can solve this equation. 

To proceed further, we need bounds on fiT as well as existence. While it is not clear that such 
bounds always hold, we can prove that they do, so long as either the Yang-Mills group G is U{1) 
(the Maxwell case), or so long as ad(V) has no globally parallel sections. We note that this last 
condition holds generically for the groups SU{2) and SU{S) of primary physical interest. The 
precise statement of the necessary boundedness result is as follows, where the Holder spaces used 
here are defined analogously to those in Definition 

Proposition 5 Suppose that either D acting on ad(V) has trivial kernel or C ~ U{1). Then 
the unique solution ^t ofD^ Dt A't = ^ Et satisfies ||Dt MtH^ q, < CT^^^ II-^'tII^. q, for some 
constant C independent of T . 
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Wc prove Proposition [SI later in this seetion. For now, let us assume the result, and proceed 
with the verification of the gluing construction conditions. To determine the size of jiT and 
therefore the size of the correction \Et — Et\, we need to estimate the quantity Dj^ Et appearing 
on the right hand side of the equation for /iy. Wc readily verify that the following bounds hold 
in the regions Q and (Recall from Section 3.4.1 that Ct\Q can be divided into the two 
components, C^^' = [-T/2,1] x S*"-! and C^^' = [1,^/2] x S"-^): 

Dt = Bi+O (e-^+''^ (46a) 
*T = *i + C'(e^*+'*) (46b) 
El = o(e-^(^)-^("-i)). (46c) 

Now, in the complement of Ct we have Et — 0, while in Q, Et — D^(£'i+i?2)- Combining 
(|46a|l . (|46b|l and H46c(l together with the identity Di *i Ei ~ wc calculate 

DtEi = *T(Di+o(e-^/2)) (*i+o(e-^/2^)Si 

A similar estimate holds for DJ, E2; thence in Q, we have Et = O (e ""-^/^) . On c!^"^ we have 
*T = *2, Dt = D2, and Et = Xi^i + i?2- Thus combining H46a|l . I)46b(l and H46c() together with 
the identity Dy *t E2 ^ we calculate 

D*tEt = *T (Di+O (6-*+")) +0 (e^)) 

= dxi A (*i + O (e~*)) El + XI *t (Di +0 (e"*+')) (*i + O (e^)) 

with an analogous estimate on c'^\ So finally, we obtain Et ^ O (e-"^/^) , from which it 
follows that 

\\ET-ET\\k,.=0(T^'\r-'^'^). 

Each of the remaining error estimates now needs to be computed. The techniques are not 
different from those used above or in the proof of Proposition |2 Hence we just summarize the 
results of the computation in the following table. 



Estimate 


S\Ct 


Ct\Q 


Q 


V)* Et 





g-nT/2 


g-nT/2 


lMT1fclM2l 


y5/2g-nT/2 


jn5/2g-nT/2 




lETlfclE2l 


y5/2g-nT/2 




e-^(^) 


lN2lfclN3l 


y5/2g-6T/2 


^_T(mi„(l,i^)) 





We note in particular that Et and Dt satisfy the momentum and energy error estimates for any 
constants k and p for which < ^ < s. and < P < min (^^, ^^) . 

The only remaining gluing conditions which need to be verified are ljNl|l and (|N4|I . Since 

ni;,D(l) = -4|i?|V(l-g)|i?D|' 
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it easily follows that for any choice of initial data, 



which establishes ljNl|l . As well, it is clear from the form of n' that IIN4|l is satisfied for any 
< c < 1. Hence we conclude that a given set of Einstein- Yang-Mills initial data can be glued at 
a chosen pair of points so long as the CK nondcgcncracy condition holds for the metric, so long 
as, if the manifold is compact, cither K ^Q,E ^Q, or Bo ^ 0, and so long as the hypotheses of 
Proposition [SI arc satisfied. 

We now return to the proof of Proposition|31 We do this first (Lemma|2l fo^' the case in which 
ad(S) has no parallel sections, and then (Lemma01) for the case in which G = U{\). The key 
step (at least when no parallel sections are present) is the establishment of lower bounds for the 
principle eigenvalue of the operator D^Dy. 

Lemma 3 Suppose ad(S]) has no parallel sections. Then there exists a constant C such that for 
T sufficiently large the lowest non-zero eigenvalue At o/D^Dt on sections o/ad(ET) satisfies 



Proof: If the statement were false, there would exist a sequence of sections Uk on Ex^. satisfying 
D^^ Dtj. Uk = AfeMfe with \k < -j^- Without loss of generality, we can assume sup|Mfe| = 1. 
Reducing to a subsequence and relabeling, we conclude using elliptic regularity that Uk converges 
uniformly on compact subsets of S* to a solution u of the equation D* Dc u = 0. 

We claim that u is not identically zero. Suppose that Uk converges uniformly to on Y,t \ Ct- 
We will show that Uk also converges uniformly to on Ct, which contradicts the assumption 
sup \uk\ = 1. First, we note that 



where := supg^^ \uk\- The volume element dV on Ct satisfies the the condition ci dVb < 
dV < C2 dVa for constants ci and C2 independent of T, where dVa is the volume element of the 
round cylinder. Letting St be the spherical cross-section of Ct at length parameter value t, we 
find that 



At> 7^. 




where the derivatives are meant in the weak sense. Hence, by Holder's inequality 





Now Kato's inequality implies that 
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and hence 

/ \uk\' < CTXk f \ukf dV + Cml 

JSt J^T 

< CTAfeVol(ST) + Cm^. 
Since VoI(St) < CT and since (by hypothesis) Afe < it follows that 

/j«r<c(l + „4), 

and we conclude from interior Schauder estimates applied to D* Dc that in fact 

sup < C ( ^ + ml 

Since converges to 0, this proves that converges uniformly to on Ct, which is the desired 
contradiction. We conclude that u is not identically 0. 

We now show that if this non zero limit u exists, then it extends to a nontrivial solution of 
D It = on S, which contradicts the hypothesis that such parallel sections do not exist. If 11 is a 
compact subset of E*, then 

|Ditp= lim / I^TfeWfel^ < liminf / iDr^Ukf = MmM Xk [ \ukf<^, 

since \uk I < 1, Vo1(StJ < CT, and Afc < j^. It follows that |D u\ = and we conclude that 
D M = on E*. Since u is bounded, it extends to a weak solution of D* D u = on E and hence 
M is a nontrivial solution of D u = on E, which is a contradiction. We therefore conclude that , 
for T large enough, Afc > □ 

The conversion of the eigenvalue estimate of Lemma |3 to the Holder estimate of Proposition 
El proceeds in a smilar fashion as was done for the vector Laplacian in |18j . From Lemma |3 
we have ll/iTlU^ < CT'^\\D^Et\\l'^- Since the volume of Et grows linearly in T we obtain 
\\D^Et\\l^- < T^/^WD^^EtWco." Finally, local Schauder estimates for D^Dc yield ||/iT||c2.° < 
C(||I?^£^t||co-° + IImtIIl^)- Hence we conclude 

\\DTtiT\y.. < CWfiriy.. < CT^mErWco.^ < CT^/^\\ET\y,.. 

The estimate of Proposition O for higher order spaces C*''" follows from another application of 
local Schauder estimates. 

LemmaOdoes not apply to U{1) bundles, which do admit parallel sections. In this case, D*D 
is the Hodge Laplacian d*d acting on scalar fields (or equivalcntly — A using our earlier notation 
and sign convention). The following Lemma implies Proposition [S] for U{1) bundles. 

Lemma 4 Any solution of d*diiT = —d*ET satisfies WdfxxWk ^ ^ ^-^"'^^^ I l^^^l Ifc q some 
constant C independent of T . 

Proof: Multiplying both sides of the equation d* d^ij^ ~ d* Et by /iy and integrating by parts we 
have 



I Idurl^ dVT= f {E,dfiT) dVr. 
The Cauchy-Schwartz inequality then implies 

IM/^tIIl2 < II-Bt||l2. 
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Since the volume of Et grows linearly we obtain 

\Wt\\l-<CT^'^\\Et\\co.^. 

This estimate replaces the eigenvalue estimate of Lemma |3 But to use effectively we have to 
show that d/iT satisfies an appropriate elliptic equation. 

Applying the Hodge Laplacian d* d -\- dd* to dfiT we have 

{d* d + dd*)dfiT — dd*dfiT ~ dd*E. 
So from local Schaudcr estimates applied to the Hodge Laplacian we obtain 

|M/iT|b.= < C{\\dd*ET\\co.^+\\dflT\\L-) 

< C{\\ET\y..+T^/^\\ET\M 

< CT'^^\\ET\y... 

The proof of the lemma in the case = 2 is complete, and the the higher order estimates follow 
from another application of local Schaudcr estimates. □ 



4.3 Einstein- Vlasov 

The Einstein- Vlasov system is used in general relativity to model self-gravitating systems of 
coUisionless matter; i.e., matter which interacts only by means of the collective gravitational 
field. There has been a renewed interest recently in establishing rigorous results for the dynamics 
of solutions of the Einstein- Vlasov system. Two useful survey papers on the subject are given in 
Em and [T|. 

On a Lorentz manifold (M, g) the additional field specified by the Einstein- Vlasov system is 
the distribution function representing the density of particles with a given space-time position 
and a given momentum. Each particle is assumed to travel along a time-like future directed 
geodesic, so its momentum at each point is mv where m is the mass of the particle and w is a 
future pointing, unit time-like vector. For simplicity we assume that all the masses are taken 
to be one; however the theory can easily be adapted to account for continuous, non-constant, 
masses. The collection 

V ~ {{x, v) \ x € Af, V a future pointing, unit timclike vector} 

forms a Riemannian hypersurfacc, called the mass shell, in the tangent bundle of M. The distri- 
bution function is then given by a non-negative function 

f -.V 

We assume for simplicity that / has compact support (again this can be relaxed). The Einstcin- 
Vlasov system is then 

where we have again chosen units so that the speed of light and Stt times the gravitational constant 
are one. The Einstein- Vlasov stress-energy tensor is given by 

Tfi„{x) ^ - f{x,v)Vf^v^dvg 

for each x € M, where dvg is the induced Riemannian volume measure on the fibre Vx of V over 

X. 
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The non-gravitational initial data for the Einstein- Vlasov system now consists of a (compactly 
supported) function /o on the tangent bundle of S. The energy density and current density of 
the non-gravitational field are then given by 

p{x) = f fo{v){l + \v\')dv, 

where 7 is the Riemannian metric on E and dv^ is the induced Riemannian volume form on the 
tangent spaces T!^^^- Note that there are no additional non-gravitational constraints which need 
to be satisfied. 

One readily checks that that if, under the conformal change 7 = "0*7, we set 

/o(«) = ^-h-^jo{<j,iv) 

then 

Ja = (f)^"^^^ Ja and p = (l>-^'^-^p. 

From this it follows that n((/)) = (j)'^^^ p = (p^i^^p and we observe that the construction now 
proceeds exactly as in the Einstein-perfect fluids case. 



5 Correcting ax in Higher Dimensions 

In section [3. 31 the following Lemma was left unproved. 

Lemma 1 Suppose there are no conformal killing fields that vanish at points pj ofT,. Then for 
T sufficiently large and for each X € C'°'"(I]t) there is a unique solution W E C'^+^'"(Sr) to 
LW = X. Moreover, there exists a constant C independent of W and T such that 

||^||fc+2,a <ct3||vf|U,„. 

The proof of the lemma is identical to the case when n = 3 found in so long as one can 
establish the following lower bound on the size of the smallest eigenvalue of the vector Laplacian. 

Theorem 2 For T sufficiently large, the lowest eigenvalue Ag = Ao(r) for L on satisfies 
Ao > CT^^ for some constant C independent of T . 

In the case n = 3, Theorem |21 follows from a perturbation argument for the lowest eigenvalue 
of the vector Laplacian on the round cylinder. Here again, the proof when n > 3 is identical to 
the case when n = 3 once one has obtained specific estimates for the vector Laplacian on the 
round cylinder, which we turn our attention to now. 

For convenience we take the vector Laplacian to operate on 1-forms rather than on vector 
fields. Let X be a covector field on the cylinder R x S"""^, which we write as 

X ^ fds + Y{s) 

where Y{s) is a covector field on S*"^^. The vector Laplacian applied to X can then be written 
in terms of its action on / and Y . We obtain 



LX^L 



f\^(^rdl + h^e ^dsSo \ff 



V ^9sde ~^di + ^6gde + ^deSe + i2-n)J \Y 
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An orthonormal basis of 1-forms on S"" ^ is given by the eigenfunctions of the Laplacian. 
Moreover, if is a 1-form such that 

A7] = Ary, (47) 
then from the Hodge decomposition and the topology of the sphere it follows that 

= fi0 + (48) 

where is a function, -0 is a divergence free covector field, and both are eigenfunctions of the 
Laplacian with eigenvalue A. So if {4>j} is an orthonormal basis of eigenfunctions of the scalar 
Laplacian with eigenvalues A,-, we have an orthonormal basis of 1-forms given by {-^d(l)j}U{ipj} 

where Sipj = and Atpj = Hjipj. Finally, we note from ^31 that for each j there exists k,l ^ N 
such that Xj = k{k + n — 2) and = (/ + 1)(Z + n — 3). 

For a covector field X of the form u{s)(f>jds + v{s)^^do4)j the vector Laplacian acts on the 
column vector (u, u)* via 



except when j = 0. In this case (po is constant, X = u{s)ds, and we simply have 

n 

Finally, for covector fields of the form X = ■w{s)ipj we have 

i-=-^5,^ + (f + 2-n). 

So together 

L = Lo0^moL'^). 
i>i 

An analysis parallel to that of ^S] shows that the temperate solutions of LX = on the 
cylinder (i.e., those with slower than exponential growth at both ends) are spanned by ds, sds, 
LUij and scoij, where for i < j, ujij = Xidxj — Xjdxi is the generator of a rotation on the sphere. 
So there is a n{n — 1) + 2 dimensional family of temperate solutions, and the bounded ones are 
also conformal Killing fields on the cylinder. 

The remainder of the differences in the analysis are contained in the following Propositions, 
which are analogous to Propositions 2 and 3 of |18| . 

Proposition 6 Let Xq — Xo{T) denote the first Dirichlet eigenvalue for L onCx — [~T/2,T/2]x 
Then 

Ao(r)>^ 

for some constant C independent of T . 

Proof: The estimate is obvious for Lq and for i" when ~ 2n{n — 2) (i.e. when i" = — ^9^). 
On the other hand when Hj = (/ + !)(/ + n — 3) for Z > 2, the lowest Dirichlet eigenvalue of L'j 
on Ct converges to ^ — (n — 2) > 0. So it remains to consider the operators L' . 



27 



Let X denote a lowest eigenfunction for L'^ with components Then 



-T/2 ^ 



2ri 



A, 



,n~l 



> 



T/2 



T/2 1^ 2 n V 2 



2 + (^^Aj + (2 - ds 



2 - n /(u')- . u 



'^2 „,2 
A 



' 2 



(50) 



2 + (^^Aj + (2 - n))v^ ds 
!i^(y')2 + i(^y'f ^ ^^2 ^ (!i^A, + (2 - n))7;2 ds, (51) 

.T/2 ^ '1 '1 

where we have integrated by parts and apphed Young's inequaUty. Since Aj > n — 1, it follows 
that + (2 - n) > i. Hence, (L^-X, X) > ^ {X, X). □ 



Proposition 7 Suppose LX ~ ^X on Ct- Then, forT sufficiently large, there exists a constant 
c independent of T such that if ^ < ^ and if Jg„_.,{\X{~T/2,0)\'^ + \X{T/2,0)\^)de < Ci, then 
for any a € [—T/2 + 1, T/2 — 1] we have 



\X{s,e)\^ dsdO < Ca 



ra+l ^ 
/a-1 JS"-2 

where C2 depends only on Ci hut not T or a, and C2 ^ as Ci 0. 



Proof: Following the proof of the corresponding result in ,1B;, it is enough to prove the result 
for each component of the separation of variables decomposition. Moreover, the estimate is clear 
for Lq and i", so we restrict our attention to L'y 



We can write L'^ in the form 



L'.: = -Adl + Bds + C 



where A, B, and C can be determined from equation H49|l . Write X = {u, v)* and /(s) ~ {AX, X) 
It follows from the equation L'^X = fiX that 

^d'jis) = \A^'^X' - i^-i/^sxp + {DX,X), 

where 



D 



_l/8(!^;^ + l/2A,-;U 

-1/16 ^r'l?^' + ^^^^^ + 2-n-^i 

I (n— l)n n ^ 

Clearly -Dii > 3Aj/8 — ji. Also, D22 can be rewritten 



(n - l)Aj 



1 fn~2 
1 - — 



16 V ?7 — 1 



+ 2 - n - ^ > 



15 [n - l)Aj 

16 n 



+ 2 - n - ^. 



When \j is 2n (its second non-zero value) we have 



15 (n- 1)A 

16 n 



771+1 



' + 2-71-^ = ^, 
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which is positive when T is sufficiently large (forcing /i to be sufficiently small). Since D22 is 
increasing in \j we have therefore proved that /(s) is convex and hence the norm of X at 
s = ±T/2 controls the Lr' norm of X over any strip a — 1 < s < a -\- 1 except possibly when 
Aj = n-l. 

To handle the case = ?i — 1 we proceed as in Proposition El above. The steps leading to 
estimate (|51|) follow as before, except now we pick up boundary terms from the integration by 
parts. Specifically, we have 

(m-fi)x,x) d9ds = / !^(u')2 + A(„')2+ 

J-T/2 n n \ n J 



H\v^ ds + b(T/2) - b(-T/2) 
n ' 



+ - - 1^1^ d^ds < b{T/2) - b(-T/2). (52) 



-T/2JS"-i J-T/2 
I -T/2 JS 

Let x(s) be a cutoff function equal to for s < and equal to 1 for s > 1, and let xt{s) 
X(s - r/2). Then 

{AX'{T/2),X'{T/2)) dO = 

^{AX\xtX\T/2)) ddds 
T_iJs^-i as 

(AX",xtX'{T/2)) + (AX',XtX'{T/2)) dOds 

— 

= / / {BX' + {C~^l)X,XTX'{T/2)) + {AX',^TX'{T/2) dOds 
< - /' / \X'\^ + \X\^ dOds + e [ \X'{T/2)\de 



where the constant c is independent of T . Taking e sufficiently small we obtain 

/ \X'{T/2)f <c [ \X'f + \xfdeds (53) 

where c is independent of T . Finally we note that 

\b{T/2)\<- [ \X{T/2f d0 + e [ \X'{T/2)f d0 (54) 

for any e > 0. Similar estimates hold at s = —T/2, and combining H52() . (|53|l . and ()54(l we 
conclude there exists a constant c independent of T such that 

I (-- 1^1' < c j \X{-T/2f + \X{T/2)t de 
J-T/2 Js"-^ \n J Js"-^ 

for /i sufficiently small, which completes the proof. □ 

The remainder of the proof of Theorem |21 now follows exactly as in ^SI and the reader is 
referred there for details. 
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6 Conclusions 



Our discussion in section 4 of the Einstein-fluid, Einstein- Yang-Mills, and Einstein- Vlasov cases 
provides a sample collection of Einstein-matter field theories for which solutions of the relevant 
constraints can generally be glued together, provided the solutions satisfy mild nondcgencracy 
hypotheses. The same is likely true for a wide collection of other Einstein-matter field theories 
as well; one need only check that the conformal method can be applied, and that the the various 
criteria discussed in sections 2 and 3 - and summarized in Theorem 1 - are met. 

It is not always easy to check these criteria. Indeed, for the Einstein-Klcin-Gordon theory, 
with the standard spacetime action principle S[g, x] — + sl^Xp + 5"^^lxP): where x is a 

C— valued scalar field, one runs into serious difficulty (even in the massless case when to = 0). 
We note, however, that for this Einstein-matter model, it is also not straightforward to use the 
standard conformal method to construct solutions of the constraint equations. The difficulty 
arises in verifying solvability for the Einstein-Klein-Gordon version of the Lichnerowicz equation 

= - 2\Dx\')i^ - ^[(^ + ^^)' + 2|^I']V'-' + i^r^ ~ 1^,^1x1')^^', (55) 

where P is a C-valued scalar field, representing the time derivative of x- For solutions of the 
constraint equations of other field theories, such as the Einstcin-Dirac theory |16| . one can readily 
check whether the conditions needed for the gluing construction presented here are satisfied. 

For those Einstein-matter theories which do satisfy our gluing criteria, can we proceed further 
and obtain stronger gluing results of the sort established in |H1-[SI ('CIP-gluing') for the vacuum 
case? To be able to show this, one needs to prove that some appropriate version of the Corvino- 
Schoen gluing results ^21 holds for the field theory of interest. Corvino appears to have done this 
for the Einstein-Maxwell theory so we should be able to obtain CIP-gluing results in this 
case. 
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